Equation of motion approach to non-adiabatic quantum charge pumping 



Amit Agarwal and Diptiman Sen 

Centre for High Energy Physics, Indian Institute of Science, Bangalore 560012, India 

(February 6, 2008) 

f*"«« , We use the equations of motion of non-interacting electrons in a one-dimensional system to nu- 

merically study different aspects of charge pumping. We study the effects of the pumping frequency, 
amplitude, band filling and finite bias on the charge pumped per cycle, and the Fourier transforms 
of the charge and energy currents in the leads. Our method works for all values of parameters, 
and gives the complete time-dependences of the current and charge at any site of the system. Our 
results agree with Floquet and adiabatic scattering theory where these are applicable, and provides 
support for a mechanism proposed elsewhere for charge pumping by a traveling potential wave. For 
non-adiabatic and strong pumping, the charge and energy currents are found to have a marked 
asymmetry between the two leads, and pumping can work even against a substantial bias. 

PACS number: 73.23.-b, 72.10.Bg, 73.63.Nm 
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I. INTRODUCTION 



The idea that periodically oscillating potentials applied at certain sites of a system can transfer a net charge per cycle 
between two leads (which are at the same chemical potential) has been studied for several years, both theoretically 
[1-18] and experimentally [19,20]. Theoretical studies have used adiabatic scattering theory [5,6], Floquet scattering 
theory [8,9] and variations of the non-equilibrium Green function (NEGF) formalism [10-12]. While Floquet scattering 
theory and the NEGF formalism work for potentials oscillating with any frequency, adiabatic theory works only for 
' low frequencies. All these methods provide expressions for the charge transferred per cycle. However, it is not easy 
to obtain from these methods the current and charge at any site as a function of time. Detailed information like 
• this may shed light on the mechanism of charge pumping, for instance, by a traveling potential wave which has been 
observed in several experiments [20]. Further, the effects of a finite bias between the leads are not easy to study 

■ analytically unless the pumping is adiabatic [21]. In this paper, we present a numerical method for obtaining all this 
information for a system of non- interacting electrons. Our method is based on solving the equation of motion (EOM) 
of the density matrix of the system [22-24]. (For time-independent Hamiltonians, this method reproduces the results 
obtained by the NEGF formalism [25]). We will see that the time-dependences of the currents and charges depend 
significantly on the amplitude and frequency of the pumping. We will show that there is a marked asymmetry in the 
Fourier transforms of the charge and heat currents in the two leads if the pumping amplitude and frequency are large. 

■ We will also show that pumping can work even if there is a substantial bias opposing it. 

As a simple model for studying charge pumping, we consider a one-dimensional system consisting of two semi- 
infinite leads a = L,R (denoting left and right) and a finite region W (wire) lying between the two. We will model all 
three regions by lattices with electrons governed by a one-channel tight-binding Hamiltonian with the same hopping 
i—i ■ amplitude —7 on all bonds, namely, 



N-l 

H = - 7 X! ( C n+l C « + 4 c n+l) , (1) 
n=l 



where N is the total number of sites. We will consider spinless electrons here; the current of spin-1/2 electrons is 
simply twice that of spinless electrons for non-interacting electrons. The dispersion of the electrons in the leads is 
Ek = —27 cos k, where k lies in the range [— it, ir]. (We are setting the Planck constant U and the lattice spacing equal 
to unity). The two leads are assumed to have the same chemical potential fi and temperature T. Time-dependent 
potentials will be applied to some sites of the wire; that part of the Hamiltonian is given by 

V(t) = X V n (t) c\c n , where V n (t) = a n cos(wi + </>„) . (2) 

n 

The sites with these potentials will be collectively called the scattering region. 
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In Sees. II and III, we will describe how Floquet scattering theory and adiabatic scattering theory respectively can 
be used to study charge pumping in the above model. In Sec. IV, we will describe the EOM method for numerically 
studying various quantities of interest. Numerical results will be presented in Sec. V. In Sec. VI, we will summarize 
our results and point out some problems for future studies. 



II. FLOQUET SCATTERING THEORY 

Briefly, Floquet scattering theory works as follows [8,9]. The incoming electrons of energy E gain or lose energy 
in quanta of uj on interacting with the scattering region. Hence, the outgoing states are characterized by energies 
E p = E + pej, where p = 0, ±1, ±2, • • •; the energies with p ^ are called the Floquet side bands. The effect of the 
scattering region can be described by a Floquet scattering matrix S a p(E p , E ), which is the amplitude for an electron 
with energy Eq entering through lead (3 to leave with energy E p through lead a. In the leads, the propagating modes 
have energies lying within the bandwidth [—2j, 27]; only these modes can contribute to charge pumping. States with 
energies lying outside the bandwidth have wave functions which decay exponentially into the leads and hence do not 
contribute to charge transfer. The wave function of an electron coming from the left lead with an energy Eo and 
wavenumber ko (with Eq = — 27cosfco) is given by 

tp{n) = e '(feon-Bot) + r P e** - *'"-***) , (3) 

p 

at a site n far to the left of the scattering region, and 

^(«) = Yl *p e^ n - E ^ , (4) 
p 

far to the right of the scattering region, where E p = —2jcosk p , and the sums over p run over values such that E p 
lies within the bandwidth of the leads. The quantities r p and t p denote reflection and transmission amplitudes in the 
different side bands; they respectively denote the elements Sll(E p ,E ) and Srl(E p ,E ) of the Floquet scattering 
matrix, where L and R denote the left and right leads. Similarly, the wave function of an electron coming from the 
right lead with an energy Eo and wavenumber fc is given by 

iP{n) = e i(-kon-E t) + f p e *(*»»- B p*) , (5) 

p 

far to the right of the scattering region, and 

V>(n) = *P e^- k " n - E ^ , (6) 
v 

far to the left of the scattering region. The reflection and transmission amplitudes are found by writing down the 
wave functions in the scattering region, and matching coefficients of terms having the same time dependence (e ±tEpt ) 
in the Schrodingcr equation at different sites. If the oscillating potentials are weak, the reflection and transmission 
amplitudes decrease rapidly as \p\ increases; at first order in the potentials, only p = ±1 contribute. The current in, 
say, the right lead is then given by 



Ir= Q r ^■[-(\ti\ 2 + \ri\ 2 ){f(Eo,i i ,T)-f(E 1 ,^,T)} 



+ — (|*-i I 2 + |r-i | 2 ) {(f(Eo,n,T) -/(£?_!, , (7) 



where f(E,/j,,T) = l/[e^ E ^/ ksT + 1] is the Fermi function, and v p = 27sin/c p is the velocity. In the limit 
lo — > 0, we have f{E±i,fi,T) — /(Eq,[i,T) = ± Lodf(E( il fi,T)/dEo, and v p /v$ — > 1. Finally, at zero temperature, 
df(Eo,n,0)/dEo = -8(Eo-ii). 



III. ADIABATIC SCATTERING THEORY 



In the limit of the pumping frequency tending to zero, the charge transport can be related to the 'frozen' scattering 
matrix S [2,3,5]. The average pumped current in this limit is found to be proportional to the frequency u, and 
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therefore the charge pumped per cycle (of time period 2tt/lu) is independent of u>. If the scattering region is connected 
to leads which are at the same chemical potential and zero temperature, the infinitesimal charge flowing from that 
region to the a-th lead is given by 

dQa = ^ (dS S\a , (8) 

where q is the electron charge, and the 'frozen' S'-matrix is evaluated at the Fermi energy Ep = — 2j cos kp = P- 
Eq. (8) can be used to relate adiabatic scattering theory to a geometric description of charge transport [3] . For weak 
oscillating potentials V n as given in Eq. (2), one finds that the charge entering lead a per cycle is given by 

AQ a = -« ]T Im\(^ r ^) i ] I dV n V m , (9) 

7T ^ L \dV n dV m /{V m }=0 iaa J 

n>m 

where the integral is done over one cycle of the oscillation. Eq. (9) is a generalization of the formula in Ref. 
[3] and can be derived from Eq. (8) by Taylor expanding the scattering matrix to first order, S(V m ) = 5(0) + 
J2 n V n (dS/dV n ){v m }=o- I n the case of two weak oscillating potentials V\(t) and V^t), Eq. (9) shows that the pumped 
charge is proportional to the area in the space of the parameters (Vi, V2). 



IV. EQUATION OF MOTION METHOD 

For a system with a finite number of sites, we study the time evolution as follows. The density matrix of the system 
evolves according to the EOM 

p{t + dt) = e - l6 ^ dt p(t) e lfl{t)dt , (10) 

where H(t) — Ho + V(t) is given in Eqs. (1-2). The current across any bond is then obtained by taking the trace of the 
appropriate current operator with p. The current operator on the bond from site n to site n + 1 and its expectation 
value at time t are given by 

jn+1/2 = ni ( c l+i c n - 4 c "+l) 1 
and j n+1/2 (t) = Tr(p(t)j n+1/2 ) = iqi [/3„ >n+ i(i) - p n+ i in (t)}. 

(11) 

The charge transferred between the right and left leads R and L between two times can be found either by integrating 
the above expression in time, or by taking the operator 

A <? = f [ E c « c « - E c « c « ] - ( 12 ) 

and computing Tr (p(t)AQ) at the two times; these methods give the same result for the charge transferred in a cycle. 
Similarly, one can compute the energy current across any site n; the corresponding operator and its expectation value 
are given by 

e« = -h 2 (cj l+ i c ™-i - 4-i c «+i) i 
and e n (t) = Tr(p(t)e n ) = - i-y 2 [p n -i,n+i(t) - p n +i,n-i(t)], 

(13) 

in a region where the on-site potential V n is zero [22] . At zero temperature, the charge current is carried by electrons 
whose average energy is given by Ep = fi; the difference between the energy current e and the charge current j 
multiplied by p gives the heat current [5,8]. We define the heat current at site n to be e„ — {p/2){j n -i/2 + Jn+i/2)- 
In all our calculations, we take the left and right leads to have Np sites each and the wire in the middle to have Nyy 
sites; the total number of sites is N — 2Np + Nw- We set the hopping amplitude 7 = 1. For calculations in which 
the chemical potential is the same in the two leads, we take the density matrix at time t = to be given by that of 
a single system governed by the Hamiltonian Hq in Eq. (1) with sites, chemical potential p and temperature T 
(which we will take to be zero). If E a and Va(") are the eigenvalues and eigenstates of the Ho (a and n label the 
states and sites respectively), the initial density matrix is given by 
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Pmn(O) = ^2 ^a(m) t/j^(n) f(E a ,fl,T) 
a 



(14) 



We then evolve the density matrix in time and compute the current and charge transferred using Eqs. (10-12). 

We will also be interested in calculations in which the chemical potentials are not the same in the two leads; suppose 
that they are given by and pr in the left and right leads respectively. In this case, we begin by setting the hopping 
amplitude to be zero on one of the bonds in the middle of the system; we then compute the density matrices pl 
and pa in the left and right parts of the system by restricting Eq. (14) to the left and right leads respectively. The 
complete density matrix p(0) is then given by the direct sum of pl and pr. We then set the hopping amplitude on 
that bond to unity, and evolve the density matrix with the Hamiltonian H(t) as usual. 

We should note that although the net charge transferred per cycle is the same across all bonds, the detailed time 
dependence of the charge transferred looks different for different bonds. (This will become particularly clear in a later 
figure). Another point to note is that the finite length of the leads (with Nl sites) implies that the system has a 
return time Tr equal to 2Nl/vf where the Fermi velocity vp — 2'ysmkp [23]; this is the time required for an electron 
to travel from the wire in the middle to the end of either of the two leads and then return to the wire. The numerical 
results can be trusted only for times which are less than Tr. Finally, there are transient effects which last for one 
or two cycles; the effects of different choices of the initial density matrix get washed out after this transient period. 
All the numerical results presented below are therefore taken from times which are larger than the transient time but 
smaller than Tr. 




V. NUMERICAL RESULTS 
A. Effect of band filling 

Figure 1 shows the charge transferred versus the Fermi wavenumber fcj? as obtained by the EOM method when 
oscillating potentials are only applied at two neighboring sites. If the amplitudes a\ and ai are much smaller than 
vf, one can use Floquet scattering theory to second order in the a n to find the charge transferred per cycle from left 
to right in the limit uj — > 0, 

a^i 01(12 ■ IM a \ cotk F 

AQ = - q — o- sin 2 - <f>i) — — . (15) 
7 2 2 

The comparison between the results obtained numerically for two different pumping frequencies and the analytical 
expression given in Eq. (15) is also shown in Fig. 1. We see that the numerical results and second order Floquet 
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scattering theory match near the middle of the band, fci? = 7r/2, but the agreement becomes poor near the band edges 
kp = and ir. This discrepancy is due to the expansion parameter a n /vF becoming large near the band edges; hence 
second order perturbation theory breaks down, and Eq. (15) is no longer valid. We can also define an adiabaticity 
parameter e = (d / v f) / {2n / lo) which is the ratio of the time taken for an electron with Fermi velocity vf to traverse 
the scattering region of length d to the time period of the oscillating potentials; the adiabatic limit corresponds to 
e — > 0. In our case, d = 1 (the oscillating potentials are at neighboring sites), while vf — 2 at kp = n/2 (half- filling). 
Even for the largest frequency of lo = tt/5 used in Fig. 1, we see that e = 1/20. This is why the two curves with finite 
lo collapse on to the results of Floquet scattering theory near fc^ = 7r/2; this may not happen if e is of order 1, i.e., if 
either word becomes larger. 

Fig. 1 shows that the charge transferred goes through an extremum as kp approaches the band edges (0 or 7r). 
This can be understood as follows. For a weak potential (ai,a 2 small), only the first Floquet side bands (p = ±1) 
contribute to the current. From the structure given in Eq. (7), one can then see that at zero temperature, the current 
gets a contribution only from electrons which lie within the energy range [Ep — lo, Ep] initially and get excited to 
the energy range [Ep 7 Ep + lo] finally. The number of states lying within these ranges starts decreasing when one 
gets very close to the band edges, namely, when Ep — lo goes below the bottom end of the band, or Ep + lo goes 
above the top end. (We recall that the band has a finite width going from —2-y to 2-y). Hence the magnitude of the 
current starts decreasing when Ep — lo falls below —2y or when Ep + lo goes above 27. Qualitatively, this is what 
one observes in Fig. 1; the charge transferred goes through an extremum when Ep = —2jkp gets within a distance 
of lo from the top end or the bottom end of the band. The extremum occurs closer to the band edge if lo is smaller 
(71-/IO instead of tt/5). 

Note that the Hamiltonian in Eqs. (1-2) is invariant under the particle-hole transformation c„ — > (— l) n <4 and 
t — > t + tt/lo, but the current operator changes sign. Hence the charge transferred is antisymmetric about kp — ir/2 
(half- filling) as we can see in Fig. 1; hence the pumped charge is exactly zero at half-filling for any value of the 
frequency lo. 



Eq. (15) can be generalized to the case in which there are oscillatory potentials at several sites as given in Eq. 
(2). In the limit lo — > 0, the use of Floquet scattering theory up to second order in the amplitudes gives the charge 
transferred per cycle from left to right to be 



Since the amplitudes a n are positive, Eq. (16) suggests that the charge transferred will be maximized if one chooses 
4>n — <Pm to be in phase with 2kp(n — m) for all pairs of sites n and m. A simple way to ensure this is to choose 
4> n = 2kpn. This is the choice of the phases cf) n made in Figs. 2-6 all of which involve systems with oscillating 
potentials applied to 8 consecutive sites, with the same amplitude a at all those sites. Hence the potential at site n 
takes the form ocos(o;i + 2kpn); this describes a potential wave traveling with a velocity u/2kp. 

Figures 2 and 3 compare the cases of weak and strong pumping. The Fourier transform of the current (computed 
at the ninth bond to the right of the scattering region) shows that for weak pumping (a = 0.2), only a small number 
of Floquet side bands contribute to the current, while for strong pumping (a = 2), a large number of Floquet side 
bands contribute. The Fourier transform of the current therefore provides a way of distinguishing between strong and 
weak pumping. 

As an application of the EOM method to recent observations of charge pumping by a traveling potential wave 
[15,20], we present in Fig. 4 the charge transferred versus time in one cycle in the strong pumping regime, as obtained 
by the EOM method for lo = 7r/100 and 7r/200, and from Eq. (8) for the adiabatic case. We see that very little charge 
is transferred in one part of the cycle, and a lot of charge is transferred in the other part; the reason for this will 
become clear below. The charge transferred per cycle is about 1.4 for lo = 7r/100, 1.8 for lo = 7r/200, and exactly 2 in 
the adiabatic case. In this model, therefore, the charge transferred increases as the pumping becomes more adiabatic. 

Figure 5 shows the density profile, Tr (p cJjC„), in and on two sides of the scattering region at eight equally spaced 
times in one cycle for the same parameters as in Fig. 4, with lo = 7r/200. We see a larger number of electrons in 
the regions where the potential has a minimum; these electrons move along with the potential minimum. The first 
six pictures in Fig. 5 show some electrons (about 1.8 in number as indicated in Fig. 4) being transported by the 
potential minimum from the right side of the scattering region to the left; during this period, very little charge is 
transferred to or from the leads. The last two pictures in Fig. 5 show these electrons being transmitted to the left 



B. Traveling potential wave 




(16) 
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FIG. 2. Fourier transform of the current (in units of q) versus frequency for oscillating potentials at 8 sites, with a = 0.2, 
uj — 7r/10 and fej? = n/4 for a system with 638 sites (Nl = 315). The left and right insets show the current and charge 
transferred versus time. 
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FIG. 3. Fourier transform of the current (in units of q) versus frequency for oscillating potentials at 8 sites, with a = 2, 
ui — 7r/10 and fcf = 7r/4 for a system with 638 sites (Nl = 315). The left and right insets show the current and charge 
transferred versus time. 

lead, while some other electrons are entering the scattering region from the right lead. These pictures illustrate the 
mechanism of charge transfer mentioned in Ref. [20]. Note that charge gets pumped in this model even though the 
'frozen' S*-matrix is almost perfectly reflecting at all times; the current would have been very small if the potential 
wave had been stationary. 



C. Asymmetry in charge and energy currents 

Although the time-averaged current must be the same at all sites due to current conservation, the fluctuations in 
the current need not be the same everywhere. This is illustrated in Figs. 6 and 7 which show the currents at the tenth 
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FIG. 4. Charge transferred (in units of q) from left to right versus time during one cycle (with period T = 2n/u>) for 
oscillating potentials at 8 sites, with o = 2, kp — n/8, and u> = 7r/100 (dash dot), 7r/200 (dashed), and tending to zero (solid 
line). 
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FIG. 5. Electron numbers (stars) and potentials at 24 sites in and on two sides of the scattering region (sites 1 to 8) at 8 
equally spaced times in one cycle, for potentials oscillating with a — 2, u> = 7r/200 and kp = 7r/8 for a system with 1618 sites 
(Nl ~ 805). The scale on the y-axis indicates both the electron number and the potential (divided by 2) at each site. 

bond on the left and tenth bond on the right respectively of the scattering region for a model with the same parameters 
as in Fig. 3. We see that the current is distinctly more noisy on the left. This is an effect of non-adiabaticity; we find 
that it disappears in the limit uj — > 0. Qualitatively, this asymmetry occurs because the potential minimum 'picks' 
up electrons from the right lead, transports them through the scattering region with a finite velocity given by u)/21cf, 
and finally 'throws' them into the left lead with that velocity; the finiteness of this velocity may be responsible for 
the additional noise on the left. 

The asymmetry can also be seen in the energy and heat currents [5,8,22] on the left and right of the scattering 
region for the same model. Figure 8 shows the outgoing energy currents ei, and en on the tenth site on the left 
and tenth site on the right respectively of the scattering region, and the product of the chemical potential [i with 
the outgoing charge currents ji and jn on the left and right, all as functions of the pumping frequency to. (Due to 
current conservation in the steady state, ji = —Jr). The differences between the two currents, namely, ej, — fij^ and 
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Current on left of scattering region 
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FIG. 6. Fourier transform of the current (in units of q) on the left of the scattering region versus frequency for oscillating 
potentials at 8 sites, with a — 2, uj = 7r/10 and Hf = 7r/4 for a system with 428 sites (Nl = 210). The left and right insets 
show the current and charge transferred versus time. 



Current on right of scattering region 
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FIG. 7. Fourier transform of the current (in units of q) on the right of the scattering region versus frequency for oscillating 
potentials at 8 sites, with a — 2, uj = 7r/10 and &f = tt/4 for a system with 428 sites (Nl = 210). The left and right insets 
show the current and charge transferred versus time. 



e R — Mr give the heat currents on the left and right; these are shown in the inset of Fig. 8. We see that the heat 
current on the left is significantly larger than on the right. We have not attempted to determine quantitatively how 
the heat current varies with lo; the dependence is known to be quadratic for small u> [5,8]. 

We should emphasize here that our model has no mechanisms (such as electron-phonon scattering) for heat dissi- 
pation in the leads. A real system will have such mechanisms, and the heat currents calculated above will eventually 
get dissipated somewhere in the leads. 
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FIG. 8. Energy currents eL and en on the left and right of the scattering region versus frequency for oscillating potentials 
at 8 sites, with a — 2 and hp = n/4. The curves marked as /xj'l and show the product of the chemical potential with the 
charge currents jh and ju on the left and right. The inset shows the heat current on the left and right versus frequency; these 
are given by cl — [ijh and en — fJ,jn respectively. 
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FIG. 9. Current (in units of q) versus bias for different frequencies (with T = 2tc/u>) for oscillating potentials at 8 sites, with 
a = 2 and ftp = 7r/4 for a system with 428 sites (Nl = 210). 

D. Effect of finite bias 

We will now consider the effect of a finite bias on the current [21]. We take the chemical potential in the left and 
right leads to be (1l = M + 9^/2 and /ir = /i — qV/2 respectively, so that there is a bias of qV between left and 
right. The calculations are done as outlined in Sec. IV. Fig. 9 shows the effect of a bias on the current for different 
frequencies for the same model as in Fig. 3. As expected, the current from left to right increases with the bias; 
however, the current continues to be negative (i.e., flows from right to left) for a range of positive values of the bias. 
For a pumping frequency of co = 27r/100 ~ 0.063, charge can get pumped against an opposing bias of almost 0.9. We 
also observe a pronounced asymmetry between positive and negative values of the bias, and the asymmetry increases 
with the pumping frequency. Finally, note that the current is very small for small bias in the adiabatic limit; this is 
because the scattering region is almost perfectly reflecting, as was mentioned earlier. In the adiabatic limit, we have 
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computed the current analytically using the Landauer-Biittiker expression averaged over time [21]. Namely, 



where < \SnL{E,t)\ 2 > t denotes the average over one time period of the 'frozen' transmission probability from left 
to right for a given energy E. 



The EOM method provides a general way of computing the charge, energy and heat currents as functions of time 
at different sites. A knowledge of the detailed space-time dependence of currents and charges is often useful. For 
instance, our results give an insight into the mechanism of charge pumping by a traveling potential wave which 
has been studied experimentally in several systems [20]. Namely, there are more electrons in a region in which the 
potential is attractive; as this region moves with time, so do those electrons. In this way, electrons are transported 
across the wire from one lead to another. We find that non-adiabatic charge pumping by a traveling potential wave 
produces an appreciable asymmetry in the Fourier transforms of the charge and heat currents. It would be interesting 
to experimentally look for an asymmetry in the Fourier transform of the current in the different systems where charge 
pumping has been demonstrated. 

The EOM approach is quite versatile and makes no assumptions about the ranges of the different parameter values; 
the potentials and pumping frequency may be small or large, and the potentials may vary with time in an arbitrary 
way, not necessarily simple harmonic or even periodic. The EOM method can also be used to study pumping at 
resonant frequencies [8,17]. 

It may be useful to compare the EOM and NEGF methods here. The EOM method requires long leads in order 
to have a large return time, so that one has a reasonable window of time to compute various quantities of interest. 
It can be extended to the case of two- and three-dimensional leads [24] ; however the calculations involve significantly 
larger systems in those cases. The NEGF formalism involves a self-energy T,(E) which takes the leads into account 
in an exact way; hence large leads do not need to be explicitly included in the numerical computations. On the other 
hand, for time-dependent problems such as charge pumping, the NEGF formalism must work with Green functions 
which would generally depend on two time arguments through the self-energy £(i, t'). In contrast to this, the EOM 
approach is local in time, and is therefore simpler to implement numerically. One can work with the NEGF approach 
in the frequency domain, but that works well only if the potentials vary harmonically in time. The EOM approach 
would work even for potentials which vary with time in an arbitrary way. 

It would be useful to extend the EOM method so as to take into account interactions between the electrons [24]. 
For small pumping frequencies and weak interactions, one can use an adiabatic Hartree-Fock approximation; for 
instance, an on-site Hubbard interaction of the form Uc^ n ^Cn^c^ \ C n ,{ can be approximated by the time-dependent 

term U[p n ^(t) & n 

c n,T c ™^]- However, this approximation will break down if the interactions are 
strong, and also if the pumping frequency is large (in which case the effective on-site interaction may not instantly 
follow the on-site densities). We would therefore require a different technique in such situations. Recently, Tomonaga- 
Luttinger liquid theory has been used to study the effects of oscillating potentials in interacting systems [26-29] . 
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